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A  Comparison  of  Non-Linear  Regression  and  Weighted  Least 
Squares  for  Predicting  Visibility  in  Germany 

by 

L.  A.  Franklin,  P.  N.  Somerville,  and  S.  J.  Bean 
University  of  Central  Florida 

1.  INTRODUCTION 

A  goal  of  the  Air  Weather  Service  is  to  be  able  to  state  the  proba¬ 
bility  that  a  weather  element  will  have  a  value  above  a  specified  threshold 
for  any  location  at  any  time.  Many  models  have  been  developed  for  such 
weather  elements  as  visibility,  ceiling,  sky  cover,  etc.,  for  locations 
where  records  exist.  Modeling  these  elements  for  locations  where  there  are 
no  records  is  a  more  difficult  task.  Some  models  have  been  developed  by 
Bean  and  Somerville  that  require  only  knowledge  of  elevation  of  the  location 
in  question  and  the  average  of  the  elevations  20  kilometers  from  the 
location  to  estimate  visibility  probabilities  for  a  specified  month  and  hour 
period  for  Germany.  These  models  were  developed  using  the  method  of  non¬ 
linear  regression.  The  method  of  evaluating  how  well  the  models  predict 
visibility  at  data-void  locations  (i.e.,  sample  reuse)  was  also  developed 
using  non-linear  regression.  This  paper  compares  the  method  of  non-linear 
regression  with  a  new  method,  "weighted  least  squares,"  in  developing  and 
evaluating  these  models. 

2.  BACKGROUND  OF  PREVIOUS  WORK  ON  VISIBILITY  IN  GERMANY 

Bean  and  Somerville  in  AFGL-TR-81-0144  "Some  Models  for  Visibility  for 
German  Stations"  showed  that  the  Weibull  distribution  was  able  to  fit 
visibility  data  for  30  stations  in  Germany.  For  any  distance  x,  the  proba¬ 
bility  F(x)  of  visibility  less  than  x  miles  is  then  given  by  the  Weibull 
formula 

F{x)  =  1  - 


where  a  different  set  of  a  and  g  values  were  derived  for  each  month  and  each 
3-hour  period. 


The  values  of  a  and  g  were  estimated  by  choosing  the  values  for  which 
the  Weibull  cumulative  distribution  most  closely  fits  the  empirical  cumula¬ 
tive  distribution.  Let  E.(x.)  be  the  empirical  probability  (step  function) 

J  * 

that  the  visibility  is  less  than  miles  for  the  station,  using  data 
from  the  RUSSWOs  (Revised  Uniform  Summary  of  Surface  Weather  Observations). 
The  values  of  x^  are  the  following  distances  in  miles:  y,  y,  4»  1, 

2 ,  y,  3,  4,  5  and  6.  Let 

sj 

Fj  (x;  ay  gj)  «  1  -  e'aj 

be  the  Weibull  distribution  at  station  j.  Then  the  values  of  a.  and  g.  that 

J  J 

are  chosen  are  those  that  minimize  the  expression 


E  E  [E-  (x.) 
j  i  J 


Fj(xiJ  «j 


(2.1) 


That  is,  the  values  of  a-  and  g.  are  those  that  minimize  the  sum  of  squares 

J  J 

of  the  difference  between  the  empirical  and  model  probabilities  over  all 
stations  and  all  distances  for  which  data  is  available.  This  deviation  is 
done  for  each  of  the  12  months  and  8  three-hour  periods  of  interest. 

The  coefficients  a,  and  g.  may  themselves  be  dependent  upon  other 

J  J 

variables.  In  fact.  Bean  and  Somerville  in  AFGL-TR-82-0335  "Some  New 
Practical  Models  for  Visibility  for  Germany  Locations  Where  No  Records 
Exist"  found  that  elevation  of  the  location  cf  interest  and  relative  eleva¬ 
tion  at  20  kilometers  around  the  location  of  interest  provided  an  improve¬ 
ment  in  the  model.  The  coefficients  can  be  written  as 
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where  EL,  and  AE,  are  transformed  and  scaled  values  for  the  elevation  at 

J  J 

station  j  and  the  average  elevation  at  20  equispaced  locations  on  a  circle 

of  radius  20  kilometers  and  centered  at  the  station.  EL  is  the  cube  of  the 

g 

elevation  in  feet,  divided  by  10  .  AE  is  the  cube  of  the  average  elevation 

g 

in  feet,  divided  by  10  .  The  constants  y0>  Yp  Y2  an<^  <50>  6i»  are 
determined  by  minimizing  expression  (2.1).  The  presence  of  all  six  con¬ 
stants  comprises  what  has  been  called  the  "variables  model"  using  non-linear 
regression.  If  we  have  only  =  Yo  and  =  60  the  model  has  been  called 
the  "constants  model"  using  non-linear  regression  and  effectively  fits  only 
one  model  a  =  6Q  and  6  =  60  for  all  stations,  ignoring  any  geographical 
features  at  the  stations.  The  "constants  model"  has  been  evaluated  and 
discussed  for  30  stations  in  Germany  in  AFGL-TR-81-0313  "Modeling  Visibility 
for  Locations  in  Germany  Where  No  Records  Exist."  The  "variables  model"  has 
been  evaluated  and  discussed  for  60  stations  in  Germany  in  AFGL-TR-82-0335. 

The  method  of  minimizing  expression  (2.1)  has  been  the  method  of 
non-linear  regression  and  has  been  discussed  in  detail  in  AFGL-TR-80-0362 
"Least  Squares  Fitting  of  Distribution  Using  Non-Linear  Regression."  The 
technique  involves  an  iterative  solution  incorporating  an  initial  estimate 
of  the  parameters.  While  the  method  has  been  extremely  successful  in 
fitting  models  and  displays  very  robust  features,  it  is  time  consuming  to 
execute  even  on  a  mid-sized  computer. 

The  method  used  to  evaluate  the  models  formulated  is  called  sample 
reuse  and  has  been  discussed  in  AFGL-TR-82-0335.  Briefly,  sample  reuse 
takes  a  single  station  and  uses  all  the  other  stations  to  build  the  model 
and  then  calculates  how  well  the  model  predicts  visibility  at  the  omitted 
station  measured  in  the  usual  root  mean  square  sense  (i.e.,  RMS).  This  is 
done  for  each  station  for  which  you  have  data  and  hence,  for  n  stations. 
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results  in  n  times  as  many  non-linear  regressions  as  needed  to  fit  the 
original  model  to  all  n  stations  at  once. 

It  was  in  the  interest  of  trying  to  find  a  new  method  of  estimating  the 
models  that  several  methods  were  investigated  in  AFGL-TR-83-0248  "A 
Comparison  of  Several  Alternatives  to  Maximum  Likelihood  for  the  Weibull 
Distribution."  In  that  simulation  study,  the  method  of  non-linear  regres¬ 
sion  seemed  more  robust  in  that  it  was  able  to  give  better  models  when  the 
data  was  contaminated  or  when  the  underlying  true  distribution  was  not  the 
form  of  the  distribution  chosen  to  model  it.  However  the  method  of  weighted 
least  squares  showed  promise  in  that  it  gave  good  results  (although  not  as 
good  as  non-linear  regression)  and  used  only  a  small  fraction  of  the  time 
that  non-linear  regression  needed.  This  method  was  initially  suggested  by 
Major  A1  Boehm,  USAF,  and  it  was  the  purpose  of  this  report  to  compare  this 
technique  with  non-linear  regression  on  actual  visibility  data  in  Germany. 

3.  METHODOLOGY  AND  RESULTS 
3.1  The  Constants  Model 

The  resulting  RMS  of  fitting  the  "constants  model"  by  both  non-linear 
regression  and  weighted  least  squares  to  the  first  group  of  30  stations  in 
Germany  is  recorded  in  Exhibit  3.1  for  all  month  and  hour  periods.  In  that 
exhibit  for  the  "constant  model,"  non-linear  regression  always  produced 
smaller  RMS  than  the  weighted  least  squares.  The  difference  was  as  small  as 
.002  for  the  period  06-08  LST  of  October  to  as  large  as  .112  for 
period  00-02  LST  of  January  where  RMS  for  weighted  least  squares  was  three 
times  the  value  for  non-linear  regression. 

In  AFGL-TR-82-0187,  "Evaluation  of  An  Observation-Based  Climatology 
Model  for  Predicting  Visibility  for  Data-Void  Locations  in  Germany,"  the 
constants  model  was  expanded  to  a  total  of  60  stations  in  Germany  and  RMS 
from  sample  reuse  based  on  non-linear  regression  was  utilized  to  measure  its 


ability  to  predict  visibility.  The  overall  results  of  Exhibits  2.3  and  2.4 
of  that  report  are  included  in  Exhibits  3.2  and  3.3  of  this  report  along 
with  the  RMS  from  sample  reuse  based  on  weighted  least  squares.  Hour 
periods  00-02  LST  and  03-05  LST  are  omitted  since  data  from  these  pre-dawn 
hours  was  frequently  missing.  For  the  stations,  the  weighted  least  squares 
method  produced  a  larger  RMS  than  non-linear  regression  at  each  station, 
usually  between  2  and  3  times  as  large.  The  average  RMS  for  non-linear 

regression  was  .098  while  for  weighted  least  squares,  it  was  .177.  Similar 
results  hold  for  the  monthly  and  hourly  RMSs. 

3.2  The  Variables  Model 

The  variables  model  was  fit  to  all  60  stations  in  Germany  and  the 
coefficients  for  the  a  and  8  terms  as  well  as  overall  fit  RMS  and  sample 
reuse  RMS  were  reported  in  AFGL-TR-82-0335  and  displayed  as  Exhibit  2.4  of 
that  report.  It  is  included  in  Exhibit  3.4  along  with  comparable  RMS 

results  of  overall  fit  and  sample  reuse  using  weighted  least  squares. 

Whether  one  compares  RMS  of  the  overall  fitted  "variables  model"  done  by 

non-linear  regression  and  by  weighted  least  squares  or  whether  one  compares 
RMS  from  sample  reuse  from  the  two  methods  the  results  are  the  same.  In 
virtually  every  instance  non-linear  regression  gives  a  smaller  RMS,  usually 
by  a  factor  of  2  to  3.  Noteworthy  exceptions  occur  at  Grosser  Falk  and 
Bad  Kreuznach  where  weighted  least  squares  does  marginally  better.  However 
both  stations  have  previously  presented  problems  of  modeling  and  have 
unusually  large  RMS  values  under  non-linear  regression. 

Exhibit  3.5  presents  the  overall  RMS  fit  of  the  variables  model  to  the 
complete  set  of  60  stations  in  Germany  using  both  non-linear  regression  (the 
top  entry  in  each  cell)  and  weighted  least  squares  (the  bottom  entry  in  each 
cell)  for  all  12  months  and  hour  periods  3  through  8.  Non-linear  regression 
in  every  case  produces  a  lower  RMS  than  weighted  least  square sometimes  as 
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little  as  10  percent  lower  (for  October  and  hour  period  9-11)  but  usually 

|  almost  50  percent  lower.  The  overall  average  RMS  for  all  months  and  time 

periods  for  non-linear  regression  was  .071  as  compared  to  .095  for  weighted 
least  squares. 

|  Exhibit  3.6  presents  the  RMS  for  the  variables  model  using  sample  reuse 

on  the  60  stations  in  Germany  using  non-linear  regression  (the  top  entry  in 
each  cell)  and  weighted  least  squares  (the  bottom  entry  in  each  cell). 

il  Again  the  RMS  achieved  by  non-linear  regression  was  always  smaller  than  the 

RMS  achieved  by  weighted  least  squares. 

4.  CONCLUSIONS 

»  The  results  of  investigating  weighted  least  squares  as  ;  ethod  to 

derive  visibility  models  is  clear  and  consistent  with  othe  esults. 
Non-linear  regression  consistently  provides  better  models  than  weighted 

M  least  squares,  usually  dramatically  better.  This  is  "better"  in  the  sense 

of  better-fitting  models  (i.e.,  overall  RMS)  and  "better"  in  the  sense  of 
better-predicting  models  (i.e.,  sample  reuse). 

||  However,  results  from  weighted  least  squares  is  consistent  with  non¬ 

linear  regression  in  the  sense  that  relatively  low  RMS  values  derived  by  NLR 
correspond  to  relatively  low  RMS  values  derived  by  WLS.  In  particular, 

*  stations  that  are  "difficult  to  fit"  or  may  contain  "bad  data"  remain  so 
under  either  technique.  Thus  it  seems  reasonable  that  WLS  can  be  used  to 
construct  many  models  of  visibility  and  explore  them  for  (relative)  goodness 

1  of  fit,  since  it  is  much  faster  to  compute  than  NLR.  Hence  having  elimi¬ 

nated  a  vast  array  of  poor  models  by  WLS,  it  would  be  possible  to  use  NLR 
upon  these  few  models  which  WLS  points  to  as  most  promising.  NLR  can  then 

*  choose  the  "best  among  the  good"  and  give  the  lowest  RMS  and  most  accurate 
values  for  the  coefficients  of  a  and  6  in  the  Weibull  Model. 
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It  is  hoped  that  such  a  two-step  method  can  be  employed  by  the  authors 

on  data  from  Norway  to  see  if  it  is  successful. 
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Hour  Period 
(LST) 


00-02 

03-05 

06-08 

9-11 

12-14 

15-17 

18-20 

21-23 

Month 

January 

.068 

.102 

.107 

.099 

.087 

.078 

.075 

.180 

.151 

.144 

.104 

.096 

.114 

.165 

February 

.082 

.093 

.078 

.064 

.057 

.053 

.072 

.090 

.178 

.135 

.081 

.067 

.060 

.057 

.088 

.164 

March 

.048 

.086 

.066 

.057 

.053 

.041 

.051 

.055 

.106 

.120 

.069 

.061 

.062 

.045 

.064 

.091 

Apri  1 

.047 

.073 

.054 

.042 

.029 

.027 

.042 

.034 

.086 

.106 

.059 

.047 

.033 

.036 

.058 

.057 

May 

.034 

.064 

.045 

.029 

.016 

.013 

.021 

.028 

.064 

.090 

.051 

.032 

.019 

.018 

.029 

.050 

June 

.053 

mmm 

am 

.035 

.019 

.016 

.026 

.027 

.082 

mm 

Km 

.042 

.024 

.020 

.042 

.048 

July 

.045 

.074 

.052 

.031 

.016 

.012 

.026 

.031 

.078 

.104 

.059 

.036 

.022 

.017 

.047 

.052 

August 

.076 

.047 

mm 

.015 

.023 

.071 

$1121 

.113 

.091 

.058 

EH 

.019 

.030 

September 

.066 

.101 

.088 

.062 

.059 

.025 

.046 

.115 

.144 

.097 

.072 

.063 

.031 

.079 

October 

.068 

.100 

.080 

.075 

.059 

.049 

.061 

.162 

.151 

.082 

.081 

.069 

.059 

.115 

.129 

November 

.062 

'.091 

.066 

.059 

.063 

.087 

.061 

.163 

.136 

.073 

.070 

.065 

.069 

.  103 

.  147 

December 

.  083 

.094 

.082 

.078 

.075 

.076 

.075 

.105 

.198 

.146 

.088 

.083 

.084 

.083 

.  100 

.192 

Exhi 

bit  3.1 

Values  of 

RMS  for 

"constants 

Model " 

for  30 

stations 

in  Germany  using 

non- 

linear  regression  (top  entry)  and  weighted  least  square  (bottom  entry). 
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WMO 

Station 

10501 

Aachen,  DL 

10224 

Bremen,  GER 

10305 

L ingen,  GER 

10427 

Kahler  Asten ,  GER 

11157 

Aigen  Ermstal ,  OS 

11448 

Plezen/Dobra,  CZ 

10900 

Bremgarten,  GER 

10948 

Oberstdorf,  GER 

10929 

Konstanz,  GER 

11120 

Innsbruck,  OS 

11150 

Salzburg,  OS 

10893 

Passau,  GER 

10908 

Feldberg,  GER 

10921 

Neuhausen,  DL 

10512 

Nurburg,  GER 

10515 

Koblenz,  DL 

10532 

Giessen,  GER 

10542 

Hersfeld,  DL 

10658 

Kissinlen,  DL 

10671 

Coburg,  GER 

10685 

Hof,  GER 

10704 

Berus,  GER 

10727 

Karlsruhe,  GER 

10742 

Ohringen,  GER 

10788 

Staubing,  GER 

10791 

Grosser  Falk,  GER 

10805 

Lahr,  GER 

10837 

Laudheim,  GER 

10953 

Kaufbeuren,  DL 

10875 

Muhldorf,  GER 

10616 

Hahn,  AB 

10610 

Bitburg,  AB 

10614 

Ramstein,  AB 

10607 

Spangdahlem,  AB 

10384 

Tempelhof,  APRT 

10755 

Ansbach,  AAF 

10544 

Fulda,  AAF 

10869 

Erding,  AS 

10765 

Feucht,  AAF 

10618 

Baumholder,  AAF 

10626 

Bad  Kreuznach,  AAF 

10971 

Bad  Tolz,  AAF 

10714 

Zweibrucken,  AB 

10633 

Wiesbaden,  AB 

10633 

Finthen,  AAF 

10763 

Furth,  AAF 

10642 

Hanau,  AAF 

Non- 

Linear 


'  at 

(North) 

Long 

(East) 

RMS 

50.78 

6.12 

.081 

53.05 

8.80 

.048 

52.52 

7.33 

.127 

51.18 

8.50 

.274 

47.53 

14.15 

.101 

49.67 

13.30 

.145 

47.90 

7.63 

.050 

47.40 

10.30 

.147 

47.68 

9.20 

.051 

47.27 

11.37 

.081 

47.80 

13.02 

.059 

48.58 

13.50 

.066 

47.87 

8.02 

.321 

47.98 

8.92 

.108 

50.33 

6.97 

.134 

50.35 

7.60 

.102 

50.57 

8.72 

.076 

50.87 

9.72 

.135 

50.20 

10.10 

.126 

50.27 

10.97 

.123 

50.32 

11.90 

.067 

49.27 

6.70 

.125 

49.02 

8.40 

.072 

49.20 

9.53 

.092 

48.82 

12.60 

.112 

49.08 

13.30 

.275 

48.37 

7.85 

.061 

48.22 

9.93 

.094 

47.87 

10.63 

.065 

48.25 

12.55 

.091 

49.95 

7.27 

.048 

49.95 

6.57 

.056 

49.43 

7.58 

.064 

49.97 

6.70 

.057 

52.47 

13.40 

.074 

49.32 

10.63 

.071 

50.53 

9.63 

.062 

48.32 

11.93 

.060 

49.38 

11.18 

.055 

49.65 

7.30 

.176 

49.87 

7.88 

.111 

47.77 

11.60 

.123 

49.22 

7.40 

.046 

50.05 

8.33 

.050 

49.97 

8.15 

.057 

49.50 

10.95 

.090 

50.17 

8.95 

.046 

(continued) 


Weighted 

L.S. 

RMS 


.169 

.172 

.168 

.178 

.202 

.159 

.164 

.198 

.170 

.198 

.181 

.157 

.217 

.146 

.108 

.191 

.194 

.180 

.179 

.171 

.140 

.160 

.178 

.174 

.185 

.161 

.157 

.159 

.166 

.164 

.147 

.172 

.184 

.177 

.192 

.174 

.171 

.177 

.171 

.199 

.189 

.200 

.169 

.165 

.160 

.187 

.171 
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Non- 

Weighted 

Linear 

L.S. 

WMO 

Station 

Lat 

Long 

RMS 

RMS 

(North) 

(East) 

10852 

Gabl ingen,  AAF 

48.45 

10.87 

.109 

.192 

10653 

Giebelstadt,  AUX  AF 

49.67 

9.88 

.082 

.186 

10687 

Grafenwohr,  AAF 

49.70 

11.95 

.073 

.169 

10734 

Heidelberg,  AAF 

49.40 

8.64 

.041 

.167 

10752 

Illesheim,  AAF 

49.47 

10.38 

.149 

.205 

10659 

Kitzingen,  AAF 

49.75 

10.20 

.066 

.189 

10763 

Nurnberg 

49.50 

11.08 

.066 

.191 

10729 

Coleman,  AAF 

49.57 

8.47 

.114 

.182 

10657 

Wertheim,  AAF 

49.77 

9.48 

.153 

.177 

10745 

Schwaebisch  Hall,  AA 

49.17 

9.78 

.097 

.179 

10712 

Sembach,  AB 

49.52 

7.87 

.037 

.169 

10862 

Siegenberg  Gunnery 

48.75 

11.80 

.177 

.219 

10738 

Echterdingen,  ARPT 

48.68 

9.22 

.057 

.183 

Overal 1 

TOM 

.177 

Exhibit  3.2 

Overall  RMS  for  Each  of  60  Stations  for  the  Constants  Model  Utilizing  Sample 
Reuse  Based  on  Non-Linear  Regression  and  Weighted  Least  Squares 
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i 

i 


S 


Hour  Period 
UST) 


06-08 

09-11 

12-14 

15-17 

18-20 

21-23 

all 

January 

.126 

.119 

.114 

.113 

.124 

.166 

.128 

.156 

.144 

.159 

.161 

.169 

.210 

.167 

February 

.117 

.099 

.099 

.111 

.152 

.115 

.147 

.137 

.159 

.167 

.175 

.217 

.167 

March 

.103 

.083 

.094 

.116 

.097 

.145 

.205 

.211 

.253 

.195 

April 

.101 

.089 

.074 

.072 

.080 

.089 

.085 

.176 

.214 

.226 

.226 

.235 

.238 

.219 

May 

.093 

.076 

.060 

.051 

.060 

.077 

.071 

.167 

.225 

.209 

.183 

.201 

.227 

.202 

June 

.094 

.073 

.054 

.045 

.056 

.070 

.067 

.171 

.212 

.184 

.159 

.182 

.183 

July 

.088 

.066 

.049 

.039 

.065 

.061 

.164 

.194 

.137 

.146 

.177 

.165 

August 

ESI 

.051 

.080 

EH 

.159 

.155 

.167 

.202 

.161 

September 

im 

.064 

.053 

.097 

.082 

HHHI 

ESI 

.163 

.161 

B 1 

.212 

.163 

October 

.100 

.097 

.082 

EBj 

mm. M 

.107 

.112 

.142 

BSfl 

November 

.112 

.108 

.104 

m 

.137 

.140 

.174 

KB 

m 

December 

.113 

.107 

.106 

mm 

m 

.143 

.136 

.153 

Ha 

•  rn  •••••: 

m 

all 

.104 

.092 

.082 

.079 

.119 

.096 

.147 

.161 

.175 

.170 

.181 

.210 

.174 

Exhibit  3.3 

RMS  Over 

All  60  Stations 

by  Month  and  Hour 

Period 

for  the 

Constants  Model 

Utilizing 

Sample  Reuse 

Based 

on  Non-Linear  Regression 

(top  entry) 

and 

Weighted  Least  Squares  (bottom  entry). 
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Regression  Squares _ 

Rel .  Sample  Overall  Sample  Overall 
Lat  Long  Elev.  Elev.  Reuse  Fit  Reuse  Fit 

Station  (North)  (East)  (ft.)  (ft.)  RMS  RMS  RMS  RMS 
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RMS  of  Sample  Re-use  and  Overall  Fit  RMS  by  Station  Over  all  Month  and  Hour  Periods  3  through  8  for  the 
Variables  Model  Based  on  Non-Linear  Regression  and  Weighted  Least  Squares 


Hour  Period 
(LST) 


Month 

06-08 

09-11 

12-14 

15-17 

18-20 

21-23 

all 

January 

.106 

.089 

.086 

.083 

.100 

.140 

.104 

.123 

.111 

.103 

.105 

.121 

.169 

.122 

February 

.074 

.061 

.060 

.079 

.083 

.087 

.080 

.081 

.102 

.155 

.102 

March 

.075 

.066 

.056 

.049 

.066 

.092 

.085 

.096 

April 

.069 

.041 

.035 

.053 

.051 

.081 

.078 

.121 

.089 

May 

.068 

.048 

.031 

.026 

.030 

.048 

.044 

.092 

.094 

.072 

.066 

.086 

.122 

.089 

June 

.070 

.048 

.031 

.026 

.032 

.095 

.093 

.088 

.087 

.086 

July 

.071 

.023 

.099 

.082 

.052 

Esfl 

.078 

August 

.082 

.053 

.036 

.028 

.035 

.053 

.098 

.067 

.063 

.043 

.062 

■ ■ 

.084 

September 

.095 

.063 

.050 

.036 

.054 

.065 

.109 

.075 

.104 

.038 

.089 

.115 

.088 

October 

.071 

.057 

.052 

.111 

.079 

.078 

.074 

.074 

.143 

.095 

November 

.076 

.064 

.115 

.083 

.087 

.086 

.105 

.143 

.101 

December 

.087 

.079 

.142 

.097 

■39 

.101 

.102 

.112 

.171 

.116 

All 

.083 

.068 

.051 

.065 

.094 

.071 

.099 

.086 

.073 

.093 

.138 

.095 

Exhibit  3.5 

RMS  by  Month  and  Hour  for  Overall  Fit  of  All  60  Germany  Stations  Using  the 
Variables  Model  by  Non-Linear  Regression  (top  entry)  and  Weighted  Least 
Squares  (bottom  entry). 


Hour  Period 
(LST) 


Month 

06-08 

09-11 

12-14 

15-17 

18-20 

21-23 

all 

January 

.118 

.111 

.099 

.115 

.150 

.117 

.130 

.119 

.113 

.114 

.131 

.176 

.131 

February 

.102 

.089 

.081 

.081 

.097 

.132 

.097 

.096 

.094 

.094 

.112 

.162 

.112 

March 

.092 

.086 

.079 

.084 

.098 

.086 

.101 

.098 

.099 

.106 

.149 

.112 

April 

.078 

.066 

.064 

.072 

.080 

.075 

.095 

.090 

.086 

.094 

.128 

.099 

May 

.085 

.040 

.038 

.052 

.060 

.108 

.124 

.094 

.078 

.104 

.136 

.107 

June 

.085 

.058 

.041 

.046 

.057 

.074 

.062 

.106 

.114 

.120 

.088 

.122 

.108 

July 

.087 

.055 

.036 

.041 

.056 

.074 

.061 

.115 

.113 

.096 

.071 

.086 

.123 

.101 

August 

.099 

.069 

.056 

.050 

.062 

.088 

.073 

.106 

.105 

.062 

.078 

.172 

.102 

September 

.105 

.067 

.060 

.097 

.083 

.118 

.096 

.138 

.085 

mm 

.126 

.112 

October 

.096 

.085 

.078 

.078 

.126 

.096 

.102 

.084 

.094 

.091 

.117 

.150 

.107 

November 

.094 

.074 

.083 

SB 

.096 

.099 

.100 

'-fas 

mm 

.111 

December 

.104 

.102 

.096 

.097 

.107 

.151 

.111 

.117 

.114 

.116 

.115 

.124 

.179 

.128 

All 

.082 

.071 

.070 

.085 

mem 

.087 

.104 

.148 

.111 

Exhibit  3.6 

RMS  by  Month  and  Hour  for  Sample  Reuse  of  All  60  German  Stations  Using  the 
Variables  Model  by  Non-Linear  Regression  (top  entry)  and  Weighted  Least 
Squares  (bottom  entry). 
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